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Manipulation of superpositions of discrete quantum states has a mathematical counterpart in 
the motion of a unit-length statevector in an iV-dimensional Hilbert space. Any such statevector 
motion can be regarded as a succession of two-dimensional rotations. But the desired statevector 
change can also be treated as a succession of reflections, the generalization of Householder transfor- 
mations. In multidimensional Hilbert space such reflection sequences offer more efficient procedures 
for statevector manipulation than do sequences of rotations. We here show how such reflections 
can be designed for a system with two degenerate levels - a generalization of the traditional two- 
state atom - that allows the construction of propagators for angular momentum states. We use 
the Morris- Shore transformation to express the propagator in terms of Morris-Shore basis states 
and Cayley-Klein parameters, which allows us to connect properties of laser pulses to Hilbert-space 
motion. Under suitable conditions on the couplings and the common detuning, the propagators 
within each set of degenerate states represent products of generalized Householder reflections, with 
orthogonal vectors. We propose physical realizations of this novel geometrical object with resonant, 
near-resonant and far-off-resonant laser pulses. We give several examples of implementations in real 
atoms or molecules. 

PACS numbers: 03.67.Lx; 32.80.Bx; 33.80.Be 
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I. INTRODUCTION 



Manipulation of discrete quantum states has long held interest, most recently for application to quantum information 
processing In the simplest realizations one deals with a nondegenerate two-state system (a "two-level atom" [2], |3( ) 
and through pulsed resonant coherent excitation produces a specified superposition of the two states, starting from a 
single state. Following the availability of laser light sources, much attention centered on producing complete excitation, 
i.e. complete population transfer from the ground state to an excited state [3j]. More recently interest has shifted to 
the production of specific superpositions of the two states, and to producing transitions between superposition states. 
In the context of quantum information a two-state system serves as a qubit; more general tasks involve multiple 
quantum states (a qunit), possibly driven by pulsed nonresonant light. 

The goal of quantum-state manipulation is to realign the statevector *&(t) from some given Hilbert-space direction 
^(ti) at an initial time U to some prescribed state ^(tf) at a final time tf, each of these statevectors being defined 
by the set of complex- valued probability amplitudes C n (t) associated with a set of quantum states \ip n ). Expressed in 
matrix form, the goal is to obtain a propagator matrix U(i/, ij) that transforms the vector C(t) of components C n (t), 

C(t f ) = U(t f ,t i )C(t i ). (1) 

As noted in the following section, there exists an extensive literature describing analytic solutions to the two-state 
system for a variety of pulsed-excitation scenarios. There exist also numerous analytic solutions to multistate systems. 

Here we consider, as a particular generalization of the two-state system, a two-level system involving two sets of 
degenerate sublevels: a less energetic set of N states \tp n ) (n = 1, 2, . . . , N) with common energy En (the lower set), 
and a more energetic set of M states \ip m ) (m = 1, 2, . . . , M) that share energy Em (the upper set). Thus we consider 
a Hilbert space of dimension A~ + M. For definiteness we here assume that A~ ^ M. We assume that interactions 
produced by laser fields can induce direct transitions between lower sublevels and upper sublevels, but not directly 
within either manifold of states (i.e. we allow electric-dipole transitions only). 

Because the statevector has unit length at all times, any such change can be regarded as a rotation, and any 
allowable motion in several dimensions can be decomposed into a succession of two-dimensional rotations, analogous 
to the three Euler angles that define an arbitrary rotation in Euclidean space. However, length-preserving changes in 
multi-dimensional space can also be produced by a succession of reflections. When expressed in matrix form these are 
generalizations of the Householder reflections used in matrix calculus [H ■ As has been shown previously, such unitary 
operations can be implemented very efficiently using coherent excitation techniques [5, 6]; as subsequently noted 0], 
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the use of quantum Householder reflections (QHR) permits efficient quantum-state engineering of transitions between 
arbitrary superpositions. 

The QHR implementation proposed earlier 0] requires a particular multistate linkage pattern, in which a set of 
N low-lying degenerate states all link, via radiative interaction, with a single upper state - a generalization of the 
tripod linkage termed an N-pod. Such a linkage pattern occurs with a lower level having angular momentum J = 1 
(three sublevels) excited to an upper level having J = (a single sublevel) but the pattern is difficult to realize for 
more than 3 lower states. It is therefore desirable to extend the QHR technique to more general linkage patterns. 

The present manuscript describes a procedure that allows implementation of the QHR when there are multiple 
states in the upper set. Such an extension makes possible the application of QHR to arbitrary angular momentum 
states, as occur with free atoms and molecules. As we show, the propagator for such situations is not a single QHR, 
but a product of QHRs, moreover with orthogonal vectors. The expression for the resulting propagator has a clear 
geometric interpretation as the effect of a succession of reflections, i.e. coupled mirrors. 

The key to this extension of the QHR is a transformation of the underlying basis states, the so-called Morris-Shore 
(MS) transformation [1, 0, This replaces the original system, with its multiple linkages between states, by a set 
of independent nondegenerate two-state systems, thereby allowing us to utilize the considerable literature of analytic 
solutions to two-state systems and produce analytic solutions for degenerate multistate systems Q. 

In this paper we present the solution of the degenerate two-level problem in a simple closed form involving sums 
of projectors of MS states. These expressions are useful for deriving analytical solutions, generalizations of known 
two-state solutions, for systems having degenerate levels. For certain conditions, wherein the transition probabilities 
between states from different sets vanish, the propagator for states of the lower set is given by a product of QHRs 
with orthogonal vectors, each vector being a bright state; a similar property applies to the upper set. We discuss 
properties of this novel geometric object, the coupled mirrors. To illustrate the procedure we develop a useful explicit 
analytic formalism for two upper states and present some examples. 



II. THE DEGENERATE TWO-LEVEL SYSTEM 



A. Two non-degenerate states 



We consider the controlled alteration of a multistate system, expressed as a redirection of the statevector, induced 
by a set of laser pulses. The simplest example of coherent excitation occurs when there are just two nondegenerate 
states, indexed 1 and 2. The time-dependent Schrodinger equation prescribes the changes in Hilbert space as (here 
and henceforth we set K = 1, thereby making no distinction between energy and frequency units) 



ijC{t) = H(t)C(i), (2) 



where the column vector C(t) has the complex- valued probability amplitudes C n (t) as elements, C(i) = [Ci(t), C2(t)] T . 
Commonly one neglects in the Hamiltonian H(i) elements that vary rapidly compared with the characteristic response 
times of the system, the so-called rotating- wave approximation (RWA) [3j . Then the Hamiltonian matrix, with suitable 
choice of energy zero-point, reads 



H(t) = 



o *n(t) 

i-fi(t)* A(t) 



(3) 



Here the slowly- varying function O(t), the Rabi frequency, quantifies the coupling between the two states. For a 
laser-driven electric-dipole transition in an atom or a molecule, Cl(t) is proportional to the transition dipole moment 
and the electric-field envelope. The detuning A(t) measures the frequency offset of the carrier laser frequency u> from 
the Bohr transition frequency loq, A = loq — lu. A time dependence in the detuning can be introduced by both u> (due 
to chirping [HI) and luq (due to Stark and Zeeman shifts). 

When the transition is exactly resonant the associated detuning vanishes, A = 0. Then the response depends only 
upon the initial conditions and upon the time- integrated Rabi frequency (the temporal pulse area), 

A{t) = f dt'\n(t')\. (4) 



At all times the probability amplitudes are expressible as trigonometric functions of A(t); when A(t) is some even- 
integer multiple of n (e.g., a 2ir pulse), the probabilities repeat their initial values. Thus the response, for resonant 
excitation, does not depend on any details of the pulse shape. Of particular use are the n pulses, which produce 
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complete population inversion between the two states, and half-7r pulses that create an equal coherent superposition 
of the two states. Resonant pulses having precise temporal area have had wide application, most notably in nuclear 
magnetic resonance [l2j and coherent atomic excitation Q. They are now a common tool in quantum information 
processing 

A variety of nonresonant (A ^ 0) pulses also lead to exact analytic expressions for the probability amplitudes Ci (t) 
and C2{t). Amongst the soluble two-state models are the models of Rabi [13| , Landau-Zener [14| . Rosen-Zener 15 1, 
Allen- Eberly @, 111, Bambini-Berman [13, Demkov-Kunike [3, Demkov [19j], Nikitin [H, and Carroll-Hioe [2l(. 
Methods for approximate solutions are also available, such as the perturbation theory, the adiabatic approximation 
[22T |. the Magnus approximation 23], and the Dykhne-Davis-Pechukas approximation [24j]. The latter, in particular, 
is a very useful tool for deriving very accurate approximations in various cases of interest, e.g., transform-limited [25j 
and chirped Gaussian puls es l 26lj . The two-state dynamics acquires interesting new features when the pulsed field 
possesses some symmetries [27j , or when it is a sequence of identical pulses (pulse train) |28[ . 

The two-state atom provides a very basic model of coherent excitation, widely used because of the relative simplicity 
with which it can be treated mathematically and fabricated experimentally. However, with this mastery has come 
interest in more general systems, involving more than two discrete quantum states. It is with a class of these that 
the present paper deals. Specifically, we describe a straightforward technique that allows efficient redirection of the 
statevector for such a system by means of laser pulses, and the use of the many exact and approximate solutions 
available for two-state systems. 



B. Multiple states as two degenerate levels 



An important extension is the linear chain of N states, each state being linked only to its nearest neighbor. 
The simplest of this is the three-state chain, involving two independent pulses, the mathematics of which has been 
extensively reported [2^] , as have lengthier chains [3(| ■ Here we consider an extension in which each of the states of 
the two-state system is replaced by a degenerate set of states that together form an energy level. Such situations occur 
commonly when one deals with angular-momentum eigenstates, as happens with the electronic structure of atoms in 
free space. For a given angular momentum J there are 2 J + 1 magnetic sublevels that, in the absence of an external 
electric or magnetic field, all have the same energy. 

We shall consider the possibility that there be no nonzero elements of the Hamiltonian matrix linking any states of 
the same energy; nonzero couplings occur only between states from the lower set and the states in the upper set. We 
make no restriction on the number of states that connect with any single state. This generalizes the usual situation 
of angular momentum states excited by electric-dipole radiation, when any state can link to no more than three other 
states (the selection rules Am = —1,0, +1). Figure [1] illustrates the linkage pattern we consider. 

We shall assume that we can control the magnitude and the phase of each interaction (e.g. each electric field 
envelope), and that every Rabi frequency has the same time dependence, which we write as 

n mn (t) = 2V mn f(t), (5) 

where V mn is a complex- valued constant and f(t) a real- valued pulse-shaped function of time bounded by unity. We 
assume that the duration of the pulse is shorter than any decoherence time, so that the dynamics is governed by the 
Schrodinger equation. Such is the situation for excitations of free atoms and molecules by picosecond or femtosecond 
pulses [llj . 
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We shall also assume that each transition has exactly the same detuning A. In the simplest situations, those 
associated with angular-momentum states, three distinct fields can be distinguished by their polarization, while 
sharing a common carrier frequency. By utilizing three independent polarization directions of a laser pulse we can 
ensure three linkages with common time dependence. More general linkage patterns, still within the model described 
here, are also possible [id ]. 



III. EXACT ANALYTICAL SOLUTION 



A. The RWA Hamiltonian 



In the RWA the Schrodinger equation provides a prescription for the time dependence of the probability amplitudes 
C n {t). Written in matrix form it is 



*-C(t) = H(t)C(i), 



(6) 



where the column vector C(t) has C n (t) as elements. We shall arrange these elements with those of the lower set first, 
followed by the upper set. Then the RWA Hamiltonian has the block structure 



H(i) 



o v/(*y 

Vt/(t) D(t) 



(7) 



Here denotes the iV-dimensional square null matrix; the zeros signify that the lower states do not interact with each 
other and that they all have the same energy, which we take as the zero-point of our energy scale. The matrix V is 
an (TV x M)-dimensional matrix whose elements V nm are the magnitudes of the couplings between the lower and the 
upper states; V* is its hermitian conjugate. Lastly, D(t) is an M-dimensional square diagonal matrix whose elements 
are all equal to the shared detuning A(i), 



D(i) = A(i)l, 



(8) 



where I is the unit matrix. The diagonal nature of D indicates the absence of interaction of the upper states amongst 
themselves. We allow the detuning A(i) to vary with time, bearing in mind use of known analytic solutions to the 
two-state model with frequency-swept detuning. It proves useful to write the matrix of interactions V as a row vector 



of iV-dimensional column vectors \V n ) = [Vim Vin> 



V 



Vn Via 
V 21 V 22 

Vni VjV2 



, Vvnj 

Via/ 
V2M 



Vnm 



[|Vi),|V 2 ),...,|V m >]. 



(9) 



B. The Morris-Shore transformation 



A significant property of the linkage pattern shown in Fig. [T] is that a transformation of Hilbert-space coordinates 
- the Morris-Shore (MS) transformation - reduces the dynamics to that of a set of independent two-state systems, 
together with decoupled states Q. When TV > M, as is the case with Fig. [1] the N = N — M additional states 
are a part of the lower-level manifold. They have no connection with excited states, and hence they cannot produce 
excitation followed by fluorescence; they are termed dark states [31], by contrast to the states that can produce 
excitation (and thence fluorescence), the bright states. The number of bright states is the lesser of N and M, in this 
case M. Figure shows the new linkage pattern in the MS basis. 

Because, by assumption, all elements of the interaction matrix share a common time dependence, the required MS 
transformation is achieved by a constant unitary transformation 



S = 



A 
B 



(10) 



Here A is a unitary TV-dimensional square matrix (AA^ = A^A = I) which transforms the lower set of states, and B is 
a unitary M-dimensional square matrix (BB^ = B^B = I) which transforms the upper set. This transformation casts 
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FIG. 2: (Color online) Linkages of Fig. [T] transformed by the Morris-Shore transformation into a set of M independent 
nondegenerate two-state systems, with couplings A n , and a set of No = N — M decoupled (dark) states. 



the dynamics into the MS basis, with new MS bright lower states {\a m )}m=ii dark lower states {\"fk)}klii an( i upper 
states {|/3m)}m=i- The transformed MS Hamiltonian has the form 



H(t) = SH(t)S t = 



_ V/(t) 
Vt/(t) D(i) 



where 



V = AVB + . 



(11) 



(12) 



Because the matrices Wt = AVV^A^ and V^V = BV^VB^ must be diagonal (possibly after removing null rows and 
rearanging the basis states), the matrices A and B are defined by the conditions that they diagonalize VV^ and V^V, 
respectively. V^V has M generally nonzero eigenvalues (n = 1, 2, . . . , M). The A^-dimensional matrix VV^ has the 
same M eigenvalues as V'V and additional Nq = N — M zero eigenvalues. From the vector form @ of V we obtain 



VV f = 



M 



{Vi\V 2 ) 
(V 2 \V 2 ) 



(VmIVx) (V m \V 2 , 



(Vi\V m ) 
(V 2 \V m ) 

(V M \V M ) 



(13a) 



(13b) 



\AV is the Gram matrix [33] for the set of vectors \V n ), (n — 1, 2, . . . , M). Thus if all these vectors are linearly 
independent then detV^V ^ and all eigenvalues of V^V are nonzero [32|]; however, this assumption of independence 
is unnecessary. 

The MS Hamiltonian (JTTJ) has the explicit form 





r 


















Ai/W 
























A 2 /(<) • 





H(t) = 




















• A M /W 












A 















A 2 /(<) • 








A •• 























A 



(14) 



The structure of H(t) shows that in the MS basis the dynamics decomposes into sets of Nq decoupled single states 
and M independent two-state systems \a n ) <-> \[3 n ) (n = 1,2,..., M), each composed of a lower state \a n ) and an 
upper state |/3 n ), and driven by the two-state RWA Hamiltonians 



X n f(t) 
Xnf(t) A 



(n = 1,2,..., M). 



(15) 
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Each of these two-state Hamiltonians has the same detuning A, but they differ in the couplings A„. Each of the 
new lower states \a n ) is an eigenstate of S/ST that corresponds to a specific eigenvalue A^ of \A/t. Similarly the new 
upper state \f5 n ) is the eigenstate of V^V, corresponding to the same eigenvalue A^. The square root of this common 
eigenvalue, A n , represents the coupling (half the Rabi frequency) in the respective independent MS two-state system 
Wn) *-* \fin)- The No zero eigenvalues of VV^ correspond to decoupled (dark) states in the lower set (we assume 
thoughout that N ^ M; therefore, dark states, if any, are in the lower set). The dark states are decoupled from the 
dynamical evolution because they are driven by one-dimensional null Hamiltonians. 



C. The propagator in the original basis 

The eigenvectors of VV^ and V^V form the transformation matrices A and B, 



They obey the completeness relations 



A = 



(7i I 

(lN \ 

("il 
(ckmI 



B = 



(fill 
(Pm\ 



(16) 



M 



N 



\a n )(a n \ +J2\lk){lk\ = I, 

n=l fe=l 
M 

£|/W.| = I- 



(17a) 
(17b) 



The propagator u n (tf, ti) for the independent MS two-state system \a n ) <-> \[3 n ) (n — 1,2,..., M) is unitary and 
therefore can be parameterized in terms of the complex- valued Cay ley-Klein parameters a n and b n , 



Un(tf,ti) 



a n 
-ble 



(18) 



where 



+ |&n| = 1 and S = h 1 A ( i ')*'- In wha t follows we express the desired control of statevector motion 
in terms of constraints on these two-state Cayley-Klein parameters. The unimportant phase factor e~ lS originates 
from the chosen representation of the Hamiltonian |7]) , which facilitates the application of the MS transformation. In 
the interaction representation, where the diagonal elements are nullified and the detunings appear in phase factors 
multiplying the couplings, the factor e~ lS disappears. (Note: This phase factor has been omitted in Eqs. (17), (18) 
and (20) of Ref. Q, where it should have been associated with a* and 6*; however, it did not affect any result there.) 

By taking into account the MS propagators (fT8"| for the two-state MS systems, the ordering of the states (fTS)) , and 
the MS Hamiltonian (TT?|) , we write the propagator of the full system in the MS basis as 





" 1 









ai 








h 















a 2 








b 2 ■ 





u = 











a-M 








b M 






-b\e- i& 























-b* 2 e- iS ■■ 








„* p — iS . 

a 2 e 















. — h* p~ iS 








• n* p~ iS 

a M e 



(19) 
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It is straightforward to show that the propagator in the original basis U = S^US reads 

U = 



where 



M 



N 



u w = 



>NM 



>MN 



U M 



^a„|a„)(a„| + ^|7 fe )(7 fe 

n=l fc=l 
M 

^ 6n|an)(/3„|, 
n=l 

M 

-e~ lS ^2b* n \p n ){a n \, 
n=l 

M 

e- i5 5>;|/3 n )(/3 n |. 



(20a) 

(20b) 
(20c) 
(20d) 
(20e) 



71=1 



The propagator Uat connects states within the lower set, Um connects states within the upper set, and Unm and 
Umn mix states from the lower and upper sets. By using the completeness relations (|17a[) and (|17bj) we find 



Uat = I + ^2(a n - l)\a n ){a n \, 
i 

M 

l + X>* -lJI&Xfl, 



Um = 



(21a) 
(21b) 



Hence the propagator Ujv does not depend on the decoupled states \jk) (k = 1, 2, ... , No). This has to be expected 
because, owing to their degeneracy, the choice of the decoupled states is not unique: any superposition of them is also 
a zero-eigenvalue eigenstate of VV^ . Because the dynamics in the original basis must not depend on such a leeway of 
choice, the full propagator U must not depend on the decoupled states at all, and this is indeed the case. 



D. Analytical solutions for degenerate levels 



Equation (|20a[) expresses the dynamics of the degenerate two-level system in terms of the dynamics of the M 
independent non-degenerate two-state systems, each with the same detuning A(t) and pulse shape f(t) of the couplings 
but with different coupling strengths A„ . Therefore, Eq. (|20a|) allows one to generalize any analytical non-degenerate 
two-state solution to a pair of degenerate levels. Recently, such a generalization of the Landau-Zener model to two 
degenerate levels has been presented 33[. This generalization displays several interesting properties, for instance, not 
all transition probabilities between degenerate states are defined for infinite time duration. Other analytical solutions 
involving two degenerate levels have been derived for five chainwise-coupled states in M or W linkage configurations 
@. We present below another interesting aspect of the solution (f2"0"|) : its geometrical nature. 



IV. QUANTUM-STATE REFLECTIONS 
A. Coupled reflections 

Of particular significance is the special case when the Cayley-Klein parameters b n are all equal to zero, 

b n = (n = 1,2,..., A/); (22) 

then all transition probabilities in the MS basis, as well as these in the original basis from the lower set to the upper 
set, vanish, i.e. Vnm = Umjv = 0. Equation (|22"]) implies that \a n \ = 1, or 



(n = l,2,...,M), 



(23) 
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for all MS two-state systems, where 4> n are arbitrary phases. 

After substituting Eqs. (|22[) and (f2"5| in Eqs. ([2U|) , the propagator in the original basis reads 



U = 



Un 
U A/ 

M 

Uat = l + ^(e^»-l)|a„)K|, 



U A / = e 



-i<5 



,1/ 



l + 5>-*»-l)|A,)<& 



(24a) 
(24b) 
(24c) 



where the propagator Uat operates in the lower set of states and Um acts in the upper set. Taking into account that 
the coupled lower states \a n ) (n = 1,2, ... ,M) are orthonormal basis vectors, i.e. (a n \ak) — S n k, as are the upper 
MS states \(3 n ) (n = 1,2, . . . , M), we rewrite Uat and Um as the products 



M 



Ujv = I| M(a n ,</) n ), 



M 



>M 



-iS 



(25a) 
(25b) 



where 

M(i/;$ = l + (e*-l)|i/)<i/|. (26) 



These individual matrices are, in fact, generalizations of the matrices used to produce Householder reflection in matrix 
computations We refer to them as generalized quantum Householder reflection (QHR) operators [H, @, 0]- The 
generalized QHR operator (|26|) is unitary, M(i/; (f))^ 1 = M(v; </>)t = M(^; — </>), and its determinant has unit magnitude, 
dct M(^; <fi) = e 1 ^. For cf> = the QHR operator reduces to the identity, M(^; 0) — I, while for <p — tt, the QHR operator 
reduces to the standard reflection, 



M(i/ ;7 r) = \-2\v)(v\. (27) 

The orthogonality of the QHR vectors \a n ) in Eq. (|25a|) automatically ensures the commutation of the QHRs, 

[M(a n , 4> n ), M(a m , <j> m )} = 0. (28) 

Therefore their ordering in the product in Eq. (|25ap is unimportant. The same argument applies to Eq. (|25bl) . 

The importance of QHR derives from the fact that any A^-dimensional unitary matrix can be decomposed into a 
set of at most N generalized QHRs [6( • The resulting ease with which pulse sequences can be designed to realize the 
QHR therefore enables one to synthesize any desired unitary transformation of a qunit state, for example a quantum 
Fourier transform [|| or transition between any two pure or mixed qunit states 



B. Special case: orthogonal interaction vectors 



In the special case when the interaction vectors \V n ) are orthogonal, {V m \V n ) = \V n \ 2 5 mn , the MS eigenvalues and 
the MS states simplify greatly. Then the matrix VW of Eq. (|13bj) becomes diagonal, and hence its eigenvalues are 
= \V n \ {n — 1,2,..., M). Moreover, the eigenstates of V^V - the MS states in the upper set - coincide with the 
original states, 

|/?„) = \<p n ), (n = l,2,...,M). (29) 
The coupled MS states in the lower set - the eigenstates of VV^ - are readily found from Eq. (113aj) , 

K) - i^T IK) = IK). (30) 
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The propagator Um in the lower set is a product of QHRs, with the normalized interaction vectors \V n ) serving as 
QHR vectors, whereas the propagator Um in the upper set is a phase gate, 

M 

U w = Y[M{V n ,<f> n ), (31a) 

n=l 

M 

Um = e- l5 ^e-^"|^„)(^„|. (31b) 

The advantage of having orthogonal interaction vectors \V n ) is that they serve as QHR vectors. The implication 
is that in order to construct a pre-selected coupled-QHR transformation (|31a|) . the required couplings are directly 
obtained from Eq. (f3T)]) . Otherwise, in the general case of non-orthogonal interaction vectors, a set of pre-selected 
QHR vectors \a n ) {n = 1,2, ... , M), defined as the eigenvectors of VV^, demand numerical derivation of the required 
couplings from Eq. (|13a,|) . 



C. Realizations 



1. Off-resonant hyperbolic-secant pulses 

The condition (123|) can be realized with the Rosen-Zener model 0, [l5| . which assumes constant detuning and 
hyperbolic-secant time dependence f(t) for the couplings, with pulse duration T, 

f{t) = sech(i/T), (32a) 
A(t) = const. (32b) 

The independent MS two-state systems share the same detuning A and pulse shape fit), but have different MS 
couplings A„. For the Rosen-Zener model the Cayley-Klein parameters a n (n — 1, 2, M) read (EL [l5l| 



- (I + izAT) 



n r(i + A„T + IiAT)r(i-A„T + IiAT)' 
where T(z) is Euler's r- function. Using the reflection formula T(l/2 + z)T(l/2 — z) = tt/ cosirz, we find 



(33) 



,2 _ 1 sin 2 (7rA„T) 
cosh 2 (±ttAT) 



k\ 2 =i- °r 2 vr:^ - (34) 



It follows that the condition \a n \ — 1 is satisfied when A„T = I (I = 0, 1, 2, . . .). The phase <f> n of a n = e 1 ^ depends 
on the detuning A, but not on the corresponding coupling A n , and for an arbitrary integer I we find Q 



l-l 



2argJJ[AT + t(2A + l)]. (35) 



fc=0 



Hence the QHR phase <p n can be produced by an appropriate choice of the detuning A. This result shows that even 
though the couplings for the MS two-state systems are not the same, the phases 4> n of the Cayley-Klein parameters 
a n coincide, <p n = <f>. This feature is unique for the sech pulse. For other non- resonant pulses, e.g. Gaussian 25j, the 
phase 4> n would depend also on the coupling and therefore will be generally different for each MS two-state system. 



2. Resonant pulses 

For exact resonance (A = 0), the Cayley-Klein parameter a n reads a n — cos(A„/2), where the pulse area is 
An = 2A„ J™ f(t)dt. When the pulse area is A n — 2(21 + l)?r (I = 0,1,2,.. .), the phase <j> n is equal to 7r; hence we 
obtain a physical realization for the standard QHR ([27)) . This result is not resticted to the sech pulse (|32ap but it is 
valid for any pulse shape with such an area When the pulse areas are multiples of 47r, the phases 4> n vanish, 

a n — 1, and the corresponding QHRs reduce to the identity. Resonant pulses therefore do not produce variable QHR 
phases <p n , which can be used as free parameters. 
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3. Far-off-resonant pulses 



Far-off-resonant pulses provide the opportunity for easy adjustment of the phase <f> n , albeit only approximately. 
Then the condition (|2"3")l is fulfilled automatically because of the smallness of the transition probabilities in each 
MS system. Specifically, if the common detuning A exceeds sufficiently much the largest MS coupling A„, then all 
transition probabilities will be negligibly small, \b n \ <C 1. By adiabatic elimination of each upper MS state one finds 

A 2 r°° 

K ~ A / f2{t)dL (36) 

For a sech pulse the integral is equal to 2T, and for a Gaussian to \fn/2. Because each MS coupling X n is generally 
different, each phase 4> n will also be different. Any desired phase <p n , or a set of such phases, can easily be produced 
by choosing the original couplings V mn , and hence the MS couplings A„, appropriately. 



V. TWO DEGENERATE UPPER STATES 



Above, we described the dynamics of the degenerate two-level quantum system in the general case when the lower 
and upper levels had arbitrary degeneracies, N and M respectively. In this section, we will illustrate these results with 
a specific example: when the upper set consists of just two degenerate states, i.e. M = 2. This case is insteresting 
because of the possible implementations in different real physical systems, several examples of which will be presented 
below. Moreover, this special case allows for an elegant analytical treatment. 



A. General case 

We retain the notation for the lower states \i/j n ) (n — 1, 2, N), and we denote the two upper states \<p') and \ip h 
The interaction matrix (O reads 



V 



V{ V" 



V V" 

V N V N 



V'),\V")], 



(37) 



where \ V') and \ V") are A^-dimensional interaction vectors comprising the couplings between the lower states and the 
corresponding upper state. The product V'V reads 



|vr (v'\v") 



(v"\V) \v'T 

With the introduction of parameters 9 and a through the definitions 

2\(V'\V")\ 

Mg(V'\V") = cr, 

we write the eigenvalues A m and the associated eigenvectors \f3 m ) {jn — 1, 2) within the upper set as 



\2 _ 
A i,2 — 

l/Ji) = 



\V'\ 2 + \V"\ 2 \V'\ 2 ~ \V"\ 2 



2 

cos( 



2cos26» ' 



e smt 
cos 9 



(38) 

(39a) 
(39b) 

(40a) 
(40b) 



The next step is to find the MS states \ot\) and \a.2} within the lower set of states. They are the eigenstates of the 
iV-dimensional matrix, 



\A/t = |V")<V"| + |V r ")<0. 



(41) 
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which correspond to the (nonzero) eigenvalues (|40ap . We construct them as superpositions of the interaction vectors 
\V') and \V"), and find after simple algebra 



K) = — (cos9\V) -e~ i(T sm9\V")) , 
Ai 

\a 2 ) = i- (e l<T sin<9|0 + cos0|V")) . 

A2 



(42a) 
(42b) 



B. The propagators 



After we have found the explicit form of the MS states (|40b[) and (|4"2"|) , we obtain the exact form of the propagators 
Un and Um, 



U 



U 



JY 



M 



M(ai,0i)M(a 2 ,^ 2 ) 

l + (e*-l)|ai)(ai| + (e* 



l)|«2)(a2| 



-iS 



L l/3i)(/?i 



! l/32)(/3 2 |] 



(43a) 
(43b) 
(43c) 



with 4>\ and 4>2 being the phases of the Cayley-Klein parameters for the MS two-state propagators (fT8|) . 

It is easy to verify that the bright states (|42ap and (|42b[) are eigenstates of the propagator Ujv with eigenvalues 
e^ 1 and e*^ 2 , respectively. Physically this means that if the qunit starts in one of these states it will end up in this 
same state, acquiring only a phase factor. This occurs because of the conditions (|22[) and the independence of the 
different MS two-state systems. If the QHR phases are equal, <j>\ = 4>2, then any superposition of \ai) and ja^) is 
also an eigenvector of \Jn- The other eigenstates of the propagator are all degenerate, with unit eigenvalue, and they 
are orthogonal to \oti) and \ot2)- F° r a qutrit (N = 3) there is only one such eigenvector (up to an unimportant 
global phase factor) and it is proportional to \ai) x \ot2). For higher-dimensional qunits, any vector in a hyperplane 
orthogonal to \a\) and \a2) is an eigenvector of Ujv- 

In the special case when the vectors |V) and \V") are orthogonal, {V'\V") = 0, the expressions simplify. Then 
6 = and a = 0, and hence 



Ai - m A 2 - \V% 

|ai) = |a 2 ) = \V") 



(44a) 
(44b) 
(44c) 



Then the propagator Uat in the lower set is a product of QHRs, in which the interaction vectors \V) and \V") serve 
as QHR vectors, while the propagator Um in the upper set is a phase gate, 



U w = 



Um = 







(45a) 
(45b) 



C. Examples 

Following are examples of linkage patterns amongst angular-momentum states, which allow application of the QHR 
theory. 



1. Two levels, J = 3/2 <-> J = 1/2 



Figure [31 shows linkage patterns possible with arbitrary polarization between the four magnetic sublevels of J = 3/2, 
the lower set, and the two of J = 1/2, the upper set. The interaction matrix has the elements (with the Clebsch-Gordan 
coefficients included), 



V 



V6 



V5v + o 

-V2V V+ 

V- -V2V 

y/3V- 



(46) 
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FIG. 3: (Color online) Linkage pattern for the four degenerate magnetic sublevels of J = 3/2, the lower set, shown coupled 
by arbitrary polarization of electric-dipole radiation to the two sublevels of J = 1/2, the upper set. The states are labeled by 
their magnetic quantum number. 




-4 -2 2 4 

Time (in units 1/7) 



FIG. 4: (Color online) Time evolution of the numerically calculated populations of the magnetic sublevels of a J — 3/2 level 
coherently coupled to a J = 1/2 level by three polarized (er + , a~ and it) pulsed laser fields, with sech shape f(t) = sech(i/T), 
V- = Vo = V+ = 8.5T -1 and detuning A = 80T" 1 . Then 6 = tt/4, a = it, AiT = 106.3, A 2 T = 38.2, <f>i = 2.65772, 
4>2 — 0.954776. The arrows on the right indicate the values derived by the QHR theory. 



where the subscripts +, — and refer to right circular (er + ), left circular (a~~) and linear (it) polarizations. 

Figure 0] shows an example of time evolution of the populations of the magnetic sublevels in the J = 3/2 level, 
starting with all population in state |-0_ 3 / 2 ) - The conditions a n = e l ^ n (n = 1, 2) are realized approximately, by using 
large detuning from the upper J — 1/2 level. In the end of the interaction, the numerically calculated populations 
are seen to approach the values predicted by the QHR theory (the arrows on the right). 



2. Two levels, J = 2 <-> J = 1 



The linkage pattern for electric-dipole couplings between sublevels of J = 2, the lower set, and J = 1, the upper 
set will, for polarization expressed as a combination of left- and right-circular polarization, appear as two uncoupled 
systems: three states form a A-linkage, while five form an M-linkage, as shown in Fig. [5] The interaction matrix is 



V = 



VEv+ o 

V- V+ 
y/W- 



(47) 



Then tan 26 = |V+F_|/[5(| V- 



\V + \ 2 )} and a = argF + — arg^_ 
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J = l 



7 = 2 



FIG. 5: (Color online) Linkage pattern for five states in an M configuration. 




FIG. 6: (Color online) Linkage pattern for three-level ladder involving a degenerate middle level. The two ends of the chain, 
magnetic sublevels with J = have the same detuning from the three intermediate sublevels of J = 1. The various linkages 
are invoked by adjusting the direction of the polarization with respect to the quantization axis. The sublevels with J = 1 form 
the lower set, while those with J = form the upper set. 



3. Three levels, J = 0^J = 1<->J = 



The MS transformation can be applied not only to a pair of degenerate levels but also to a ladder of degenerate 
levels, as long as there is only a single detuning. Figure [5] illustrates an example in an angular momentum basis 
between the magnetic sublevels of three levels with angular momenta J = 0, 1, 0. The interaction matrix has the form 



V 



VL V". 



(48) 



and the formalism of this section applies. In this case all fields can be changed independently. 



VI. DISCUSSION AND CONCLUSIONS 



We have here extended the earlier work on QHR in Hilbert space to allow more general linkage patterns between 
the quantum states, with particular attention to degenerate sublevels that occur with angular momentum states. The 
extension relies on the use of the Morris-Shore transformation to reduce the original multi-linkage Hamiltonian to a 
set of independent two-state systems, thereby allowing the utilization of various known two-state analytic solutions. 
Such solutions, when expressed in terms of Cayley-Klein parameters, readily lead to conditions upon the pulse areas 
and the time-varying detunings of the excitation pulses. 

We have found that the propagator within the lower (or upper) set of degenerate states represents coupled quantum- 
state mirrors. The realization of these coupled QHRs within the lower (or upper) set of states requires certain 
conditions on the interaction parameters, specifically that all transition probabiities between the lower and upper sets 
must vanish. We have proposed three physical realizations of this condition: with resonant, near-resonant (hyperbolic- 
secant), and far-off-resonant pulses. Resonant and near-resonant pulses provide exactly zero transition probabilities, 
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but offer less flexibility in the QHR phases 4> n ; moreover they require a carefully chosen pulse area in each MS 
two-state system, which leads to a number of conditions on the interaction parameters. Far-off-resonant pulses fulfill 
the zero-probability conditions only approximately, but offer much more flexibility. Then the only restriction is for 
sufficiently large detuning, without specific constraints on pulse areas, because the zero-probability conditions are 
fulfilled simultaneously in all MS two-state systems. 

For angular-momentum states, there are six independent interaction parameters: three polarization amplitudes, 
two relative phases between different polarizations, and the common detuning. Therefore the constructed QHR has 
six free parameters (with the far-off-resonance realization). For a more general linkage, the number of independent 
parameters can be much larger. 

In a subsequent paper [34j , we shall describe the mathematical aspects of this largely unknown operator of coupled 
reflections, which however, as we have shown here, arises naturally in quantum systems. In particular, we shall show 
how one can factorize, and therefore synthesize, an arbitrary U(JV) propagator by such objects. The procedure is 
more efficient than a set of rotations that would produce the same transformation. 

We conclude by pointing out that the confinement of the statevector evolution to the lower set of states, and the 
availability of simple and efficient tools for its engineering, such as coupled QHRs, can be an essential ingredient for 
decoherence-free quantum computing [351 ] - 
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